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The elliptic flow in collisions of neutron-rich heavy-ion systems at intermediate energies 
emerges as an observable sensitive to the strength of the symmetry energy at supra- 
saturation densities. First results obtained by comparing ratios or differences of neutron 
and hydrogen flows with predictions of transport models favor a moderately soft to linear 
density dependence, consistent with ab-initio nuclear-matter theories. Comprehensive 
data sets of high accuracy can be expected to improve our knowledge of the equation of 
state of asymmetric nuclear matter. 



1. Introduction 

Heavy-ion reactions at sufficiently high energy represent the only means for com- 
pressing nuclear matter in laboratory experiments and for studying the nuclear 
equation of state (EoS) at supra-saturation densities.^ From the extensive search 
for observables suitable for probing the brief high-density phase of the collision, col- 
lective flows and sub-threshold production of strange mesons have appeared as most 
useful. The present consensus that a soft EoS, corresponding to a compressibility 
K w 200 MeV and including momentum dependent interactions, best describes the 
high-density behavior of symmetric nuclear matter is base d on s tudies of flow and 
kaon production within the framework of transport theoryP^El 

In recent years, motivated by the impressive progress made in observing proper- 
ties of neutron stars and in understanding details of supernova explosion scen arios, 
the EoS of neutron-rich asymmetric matter has received increasing attention.SEE] 
The symmetry energy, which is the difference between the energy per nucleon of 
neutron matter and of symmetric matter, is considered to be one of the biggest un- 
knowns in this context. We have precise information on the symmetry energy near 
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saturation density from the knowledge of nuclear masses. For densities below satu- 
ration, investigations of heavy ion reactions in the Fermi energy regime have con- 
strained the symmetry energy considerably^ and the importance of clus tering for 
the symmetry energy at very low density has recently been demonstrated EMHl At 
super-saturation density, however, the symmetry energy is still largely unknown for 
several reasons. Phenomenological forces, even though well constrained near satura- 
tion, yield largely diverging results if they are extrapolated to higher densities PEH 
Many-body calculations with realistic potentials face the problem that three-body 
forces and short-range correlations are not suffic iently constrained at higher densi- 
ties at which their importance increases fEUHHl Even the magnitude of the kinetic 
contribution, related to the nuclear Fermi motion and considered as principally 
understood, is possibly modified by a redistribution of nucleon momenta due to 
short-range correlations in high-density nuclear matterffSl 




Fig. 1. Isospin asymmetry <5 = (p n — p p )/p as a function of the normalized density p/po at time 
t = 20 fm/c in the 132 Sn + 124 Sn reaction with a stiff (E° ym ) and with a soft (E^ ym ) density 
dependence of the nuclear symmetry energy. The corresponding correla tion for neutron stars in 
/3-equilibrium is shown in the inset (reprinted with permission from Ref.ES Copyright (2002) by 
the American Physical Society). 



The appearance of the symmetry energy in many aspects of nuclear structure 
and reactions would seem to imply that constraints of various kinds should be avail- 
able. However, quantities as, e.g., the thickness of the neutron skin in heavy nuclei 
or the isospin transport in reactions of isospin-asymmetric nuclei in the Fermi- 
energy domain are predominantly sensit ive to the s trength of the symmetry term 
at densities below the saturation valueP ^^^^ Extrapolations of parameteri- 
zations to higher density raise the question of up to where they remain realistic. 
This emphasizes the need for more direct high-density probes. As in the case of 
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symmetric nuclear matter, collective flows and sub-threshold particle production 
are obvious candidates. 

A strong motivation for exploring the information contained in isotopic flows was 
provided by Bao-An Li when he pointed to the parallels in the density-dependent 
isotopic compositions of neutron stars and of the transient systems formed in colli- 
sions of neutron-rich nuclei as a function of the EoS input used in the calculation!^ 
Figure [T] illustrates the remarkable fact that the same physical laws are confidently 
applied to objects differing by 18 orders of magnitude in linear scale or 55 orders of 
magnitude in mass. It will allow us to infer properties of these exotic astrophysical 
objects from data obtained in laboratory experiments. The main difficulty is the 
comparatively small asymmetry of available nuclei. Symmetry effects are, there- 
fore, always small relative to the dominating isoscalar forces which one hopes will 
cancel in differences or ratios of observables between isotopic partners. The observ- 
able proposed by Li is the so-called differential directed flow which is the difference 
of the multiplicity-weighted directed flows of neutrons and protons. Directed flow 
describes the rapidity dependence of the mean in-plane transverse momenta of ob- 
served reaction products. 

A further encouragement was provided by transport model calculations accord- 
ing to which the elliptic flow of free neutrons and protons responds significantly to 
variations of the parameterization of the symmetry energyP^ Elliptic flow relates to 
the azimuthal anisotropy of particle emissions. This has motivated a reanalysis of 
the FOPI/LAND data for 197 Au + 197 Au collisions at 400 MeV/nucleon, collected 
many years ago a nd use d to demonstrate the existence of neutron squeeze-out in 
this energy regime Squeeze-out refers to a dominant out-of-plane emission of 
particles, relative to in-plane emission, and is considered as evidence for the pres- 
sure buildup in the collision zone. The analysis favors a moderately soft to linear 
density dependence of the symmetry energyPH 

This finding had a particular significance, in spite of a large statistical uncer- 
tainty. Rather different conclusions, suggesting either a super-stiff or super-soft 
behavior of the symmetry energy, had previously been reached in analyses of the 
7r~ /ir + production ratios, measured by the FOPI Colla boration for the same 197 Au 
+ 197 Au reactionpSI with different transport models PAUSES l n particular, the 
super-soft result has initiated a broad discussion of how it might be reco nciled 
with other observations as, e.g., observed properties of neutron stars j2fi| 28 | 29| rj ne 
FOPI/LAND elliptic- flow data were found to be inconsistent with this extreme as- 
sumption also when compared to recent QMD transport-model calculations, fairly 
independent of particular choices made for other model parameters!^ 

The obvious need to improve the statistical accuracy beyond that of the existing 
data set has initiated a dedicated measurement of collective flows in collisions of 
197 Au + 197 Au as well as of 96 Zr + 96 Zr and 96 Ru + 96 Ru which has been carried 
out in 2011 with the LANDAU detector coupled to a subset of the 

chimera^ 

detector arrayP^ I n the following sections, the present situation will be described 
in more detail. 



March 3, 2013 21:38 WSPC/INSTRUCTION FILE elliptic'ijmpc 



4 W. Trautmann and H.H. Walter 




Fig. 2. EoS in nuclear matter and neutron matter. BHF/DBHF and variational calculations with 
realistic forces are compared to phenomenological density functionals NL3 and DD-TW and to 
ChPT. The left panel zooms the low density range (from reprinted with kind permission 

from Springer Science+Business Media). 

2. Present Knowledge of the Symmetry Energy 

Our knowledge of the symmetry energy is originally based on nuclear masses whose 
dependence on the isotopic composition is reflected by the symmetry term in the 
Bethe-Weizsacker mass formula. A density dependence is already indicated by the 
use of separate bulk and surface terms in refined mass-formulae. In the Fermi-gas 
model, it is given by a proportionality to (p/po) 7 with an exponent 7 = 2/3, where 
Po ~ 0.16 nucleons/fm is the saturation density. This so-called kinetic contribution 
to the symmetry energy, however, amounts to only about 1/3 of the symmetry term 
of approximately 30 MeV for nuclear matter at saturation. The major contribution 
is given by the potential term reflecting properties of the nuclear forces. 

Nuclear many-body theor y has pr esented us with a variety of predictions for 
the nuclear equation of state ! I The examples shown in Fig. [2] for the two 

cases of symmetric nuclear matter and of pure neutron matter demonstrate that, 
overall, the results are quite compatible among each other, except for densities 
exceeding saturation at which the predictions diverge. The symmetry energy E syia 
can be defined as the coefficient of the quadratic term in an expansion of the energy 
per particle in the asymmetry S = (p n — p p )/p, where p„,p p , and p represent the 
neutron, proton, and total densities, respectively, 

E/A(p, 6) = E/A(p, 8 = 0) + E sym (p) ■ S 2 + 0(<5 4 ). (1) 

In the quadratic approximation, the symmetry energy is the difference between 
the energies of symmetric matter (S = 0) and neutron matter (5 = 1). Also the 
symmetry energy diverges at high density, as expected from Fig. [2j while most 
empirical models coincide near or slightly below saturation, the density range at 
which constraints from finite nuclei are valid (Fig. [3]). 
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Fig. 3. Symmetry energy as a function of density as predicted by different models. The left panel 
zooms the low density range up to saturat ion. The full lines represent the DBHF and variational 
approaches using realistic forces (from RefU^, reprinted with kind permission from Springer Sci- 
ence+Business Media). 



The behavior of the symmetry energy at very low densities is not correctly 
described in homogeneous mean-field approaches which do not include clustering 
effects. This is a topic of considerable current interest with go od agreement being 
reached by experimental and theoretical investigations jSEUQI The uncertainty at 
supra-saturation density in phenomenological models may be considered as the 
expected consequence of extrapolations leading beyond the range at which the 
chosen parameterizations are effectively tested!^ 

In calculations using realistic forces fitted to two- and three-nucleon data, the 
uncertainty is mainly related to the short-range behavio r of the n ucleon-nucleon 
force and, in particular, to the three-body and tensor forces .H 5 | 14 | 15 | ^ ne three-body 
force has been shown to ma ke an essential but quantitatively small contribution to 
the masses of light nucleiP^The extrapolation of the partly phenomenological terms 
used there to higher densities is, however, highly uncertain!^ The general effect 
of including three-body forces in the calculations is a stiffening of the symmetry 
energy with increasing densityPSEll Short-range correlations become also increas- 
ingly important at higher densities; results from very recent new experiment s will, 
therefore, have a strong impact on predictions for high-density nuclear matter P^ I ^ I 
Increasingl y pr ecise data from neutron-star observations may further provide useful 

loo] 

constraints. 11221 

At hig her energ ies, the momentum dependence of the nuclear forces becomes 
import ant j 5 | 6 | 39 | 40 | j^. j g we \\ k n0 wn that nuclear mean fields are momentum depen- 
dent, as seen, e.g., in the energy dependence of the nuclear optical potential. The 
dominating effect is in the isoscalar sector but there is also an important isovector 
momentum dependence. It manifests itself as an energy dependence of the isospin- 
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dependent part of the optical potential but can also be exp ressed in terms of a 

II Ol 

difference of the effective masses of protons and neutrons. 1 ^. Even the ordering of 
these effective masses is still an open problem. It has, moreover, been shown that 
the effective mass differences and the asymmetry dependence of the EoS are both 
influencing particle yields and flow observabl es. and t hat additional observables will 
be needed to resolve the resulting ambiguity! I I I 




Fig. 4. Parameterizations of the nuclear symmetry energy as us ed i n transport codes: three pa- 
rameterizations of the potential term used in the UrQMD (Ref.ED) with power law coefficients 
7 = 0.5, 1.0, and 1.5 (lines with symbols as indicated), the res ult with 7 = 0.69 obtained from 
analyzing isospin diffusion data with the IBUU04 (full line, Ref.^21), an d the super-soft and stiff 
parameterizations obtained from analyzing the 7r — Arj" production ratios with the IBU U04 (dot- 
ted line, Refl^Sl) a nd the LQMD (dashed line, Refl^ZI) transport models (from Ref.H31 reprinted 
with kind permission from Springer Science+Business Media). 



Transport theories needed for calculating the temporal evolution of nuclear reac- 
tions often use simplified parameterized descriptions of the composition-dep end ent 
part of the nuclear mean field. In the UrQMD of the group of Li and Bleicherj^Hthe 
potential part of the symmetry energy is defined with two parameters, the value at 
saturation density, usually taken as 22 MeV in their calculations, and the power-law 
coefficient 7 describing the dependence on density, 

E sym = ££ y °t + = 22 MeV ■ (p/p ) 7 + 12 MeV ■ (p/ P o) 2/3 - (2) 

In other codes the nuclear potential of Das et al. with explicit momentum depen- 
dence in the isovector secto r is used321 There, as in the IBUU04 developed by the 
groups of Li and Chenj^ESl the density dependence of the symmetry energy is char- 
acterized by a parameter x appearing in the potential expressions. Examples of 
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these parameterizations and of results obtained from the analysis of experimental 
reaction data are given in Fig. |4j The stiff (E£ ) and soft (E^ ym ) density depen- 
dences of Fig. [T] correspond approximately to the cases 7 = 1 and x — 1 shown 
here. 

Parameterizations of this kind have the consequence that, once the symmetry 
energy at the saturation point is fixed, a single value at a different density or, 
alternatively, the slope or curvature at any density will completely determine the 
parameterization. Measurements of a variety of observables in nuclear structure and 
reactions have, therefore, been used to obtain in this way results for the density 
dependence of the symmetry energy. They are often expressed in the form of the 
parameter L which is proportional to the slope at saturation, 



Most results with their errors fall into the interval 20 MeV < L < 100 MeV and are 
compatible with a most pro bable valu e L w 60 MeV, roughly corresponding to a 
power-law coefficient 7 = 0.6.^^^^ The f u vj ]j ne m Fig.|4]represents, e.g., the re- 
sult deduced by Li and Chen from the MSU isospin-diffusion data and the neutron- 
skin thickness in 208 Pb.E31 The corresponding slope parameter is L = 65 MeV. 
Rather similar constraints have been d educed fr om very recent investigations and 
observations of neutron-star properties JHEZEU 

Considerable progress regarding the correlation of the symmetry energy with 
particular observ ables in different models has been made by the Florida and 
Barcelona groups In continuation of the work of Typel and Brown and of 
FurnstahlpSGETa universal correlation between the thickness of the neutron-skin of 
208 Pb and the L parameter has been found for empirical mean-field interactions!^ 
The determination of the neutron-skin thickness by measuring the parity-violating 
contribution to electron scattering at high energy will thus offer a practically model- 
free access to the slope at saturation, even though it is obtained by probing nuclear 
matter at an average density of typically 2/3 of this valueP22 

3. Dynamics of Heavy-Ion Reactions 

Densities of two to three times the saturation density may be reached on time 
scales of « 10 - 20 fm/c in the central zone of heavy-ion collisions at relativistic 
energies of up to « 1 GeV/nucleon!^ The resulting pressure produces a collective 
outward motion of the compressed material whose strength will be influenced by 
the symmetry energy in asymmetric systems.^ At the same time, the excitation 
of A resonances in hard nucleon-nucleon scatterings leads to the production and 
subsequent emission of charged and neutral it and K mesons. The relative intensities 
of isovectors pairs of mesons depend directly on the proton-neutron content of the 

I r V)| 

site where they are produced}^ suggesting them as sensitive probes for the high- 
density symmetry energy. 

Measurements of K + /K° production ratios have been performed by the FOPI 
Collaboration for the neutron- rich and neutron-poor A = 96 systems 96 Zr + 96 Zr 



L = 3po ■ dE sym /dp\p . 



(3) 
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and 96 Ru + 96 Ru!S2l The experimental acceptances for these particles were consid- 
erably different, because of their different decay channels, and the precision required 
for an isotopic study was more reliably achieved with double ratios obtained from 
the data for the two systems with different isotopic compositions. A significant sen- 
sitivity to the symmetry energy was expected from calculations under equilibrium 
conditions but was found to be considerably reduced when actual collisions were 
modeled.^ 

The FOPI Collaboration has also presented an extensive account of the ir~ /ir + 
production ratios meas ured for the reaction 197 Au + 197 Au at several energies up 
to 1.5 GeV/nucleon|231 Theoretical analyses of these data, however, come to rather 
conflicting conclusions, sugges ting ever ything from a rather stiff to a super-soft be- 
havior of the symmetry energy!2£12212ZlThe reason may lie in the treatment of the A 
dynamics and in competin g eff ects of the mean fields and the A thresholds which 
will require further studies!^ In the analysis of Xiao et et al. with the IBUU04 
transport modelj^l the measured system and impact-parameter dependences have 
been reproduced, however, only by assuming a super-soft dependence of the sym- 
metry energy on density, close to the x = +1 case shown in Fig. [31 The presently 
inconclusive situation for pion ratios is very unfortunate because variations of up 
to 20% for soft versus stiff parameterizations are expected. 

Dynamical flow observables have been pr oposed by several groups as probes for 
the equation of state at high densityP^EH among them the so-called differential 
neutron-proton flow which is the difference of the parameters describing the collec- 

|Of|| 

tive motion of free neu trons and protons weighted by their numbers. As pointed 
out by Yong et aZ.j^this observable minimizes the influence of the isoscalar part 
in the EoS while maximizing that of the symmetry term. Its proportionality to 
the particle multiplicities, however, makes its determination very dependent on the 
experimental efficiencies of particle detection and identification and on the precise 
procedure for distinguishing free and bound nucleons in calculations. Therefore, 
also flow differences or ratios have been considered. 

4. Elliptic Flow 

The so-called squeeze-out of nuclear matter from the compressed interaction zone 
has first been observed in experiments at the BevalacP^In a sphericity analysis, the 
event shape in three dimensions was characterized by a kinetic-energy flow tensor 
whose main orientation with respect to the beam direction represents the collective 
sidewards flow and whose cross section, if it is non-isotropic, indicates the existence 
of elliptic flow. At the bombarding energies of up to 1.05 GeV/nucleon investigated 
in these studies, a preferential emission of charged particles perpendicular to the 
reaction plane has been observed. The shadowing by the spectator remnants as they 
pass each other during the collision reduces the in-plane flow, so that the strength 
of the off-plane emission, as quantified by the azimuthal anisotropy, reflects the 
internal pressure. 
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Fig. 5. Elliptic flow parameter V2 at mid-rapidity for 197 Au+ 197 Au collisions at intermediate 
impact parameters (about 5.5-7.5 fm) as a function of incident energy, in the beam frame. The 
filled and open circles represent the INDRA and FOP I datajSSES! respectively, for Z = 1 particles, 
the trian gles represent the EOS and E895 datJSZ] for protons and the square represents the 
E877 dataSSI for all charged particles (from 

RefEU 

reprinted with kind permission from Springer 

Science+Business Media). 



It has become customary to express both, directed and elliptic flows, and pos- 
sibly also higher flow components by means of a Fourier decomposition of the 
azimuthal distributions measured with respect to the orientation of the reaction 

plane ^EEHSSI 

dN N ( v-^ ,| 

d(^M = ^ y 1 + 2 £ w ™ cos - **> J • (4) 

where iVo is the azimuthally integrated yield. The coefficients v n = (cosn(</> — 4>r)) 
are functions of particle type, impact parameter, rapidity y, and the transverse 
momentum pt; V\ and V2 represent the directed and elliptic flows, respectively. 

Elliptic flow has become an important observable at other energy regimes as 
well. At ultrarelativistic energies, the observation of the constituent-quark scaling 
of elliptic flow is one of the prime arguments for deconfinement during the early 
collision phase, and properties of the formed q uark-gluon l iquid are deduced from 
the observed magnitude of collective motions ISQIfill£2lfiS I it implies that elliptic 
flow develops very early in the collision which is valid also in the present range of 
relativistic energies, as confirmed by calculations P^Isotopic flow differences appear 
thus very suitable for studying mean-field effects at high density. 

An excitation function of the elliptic flow of Z = 1 particles in 197 Au+ 197 Au 
collisions from various experiments is shown in Fig.[5j Squeeze-out perpendicular to 
the reaction plane, i.e. 1)2 < 0, as a result of shadowing by the spectator remnants is 
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observed at incident energies between about 150 MeV/nucleon and 4 GeV/nucleon 
with a maximum near 400 MeV/nucleon. At lower energies, the collective angular 
momentum in the mean-field dominated dynamics cause the observed in-plane en- 
hancement of emitted reaction products. The figure also shows that elliptic flow 
can be measured quite precisely. The reliability of the applied methods is being 
demonstrated by the good agreement of data sets from di fferent exp eriments in the 
overlap regions of the studied intervals in collision energyP^EHZOT 
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Fig. 6. Measured flow parameters jji (top) and V2 (bottom) for mid-peripheral (5.5 < b < 7.5 fm) 
197 Au + 197 Au collisions at 400 MeV/nucleon for neutrons (dots), protons (circles), and hydrogen 
isotopes (Z = 1, open triangles) integrated within 0.3 < pt/A < 1.3 GeV/c/nucleon as a function 
of the normalized rapidity yi^h/l/p- The UrQMD predictions for neutrons are shown for the FP1 
parameterization of the in-medium cross sections and for a stiff (7 = 1.5, full lines) and a soft 
(7 = 0.5, dashed) density dependence of the symmetry term. The e xper imental data have been 
corrected for the dispersion of the reaction plane (reprinted from Refl^AJ Copyright (2011), with 
permission from Elsevier). 




5. Results from FOPI/LAND 

The squeeze-out of neutrons has first been observed by the FOPI/LAND Collabora- 
tion who studied the reaction 197 Au + 197 Au at 400 MeV/nucleon!^ The squeeze- 
out of charged particles reaches its maximum at this energy (Fig.[5|), and similarly 
large anisotropies were observed for neutrons P^l The neutrons had been detected 
with the Large- Area Neutron Detector, LANDAU while the FOPI Forward Wall, 
covering the forward hemisphere of laboratory angles 6*i a b < 30° with more than 
700 plastic scintillator elements, was used to determine the modulus and azimuthal 
orientation of the impact parameter. 
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The motivation for returning to the existing data set has been provided by stud- 
ies performed with the UrQMD transport code for this fairly neutron-rich system 
(N/Z = 1.49) which indicated a significant sensitivity of the elliptic-flow parame- 
ters to the assumptions made for the density dependence of the symmetry energy!^ 
In calculations with power-law coefficients 7 = 0.5 and 1.5 (cf. Eq. [2] and Fig. [4]), 
the relative strengths of neutron and proton elliptic flows were found to vary on 
the level of 15%. 

The performed reanalysis of the data consisted mainly in choosing equal accep- 
tances for neutrons and hydrogen isotopes with regard to particle energy, rapidity 
and transverse momentum (energy and momentum per nucleon for deuterons and 
tritons). The results obtained for a mid-peripheral event class are shown in Fig. [5] 
as a function of the rapidity y, normalized with respect to the projectile rapidity 
y p = 0.896. Their asymmetry with respect to mid-rapidity, y/y p = 0.5, is caused by 
the LAND detector whose kinematic acceptance in p t increases with y. With the 
yields dropping at large pt , the statistical errors become large at forward rapidity. 
The theoretical predictions have been obtained simulating the LAND acceptance 
and the experimental analysis conditions. The results, shown for neutrons in Fig.Hl 
follow qualitatively the experimental data. 

According to the UrQMD model, the sensitivity of the directed flow of neutrons 
to the stiffness of the symmetry energy is nearly negligible (Fig. [6l top panel) 
while it is significant for the elliptic flow (bottom panel). This is evident also from 
the dependence of the elliptic flow parameter vi on the transverse momentum per 
nucleon, pt/A, shown in Fig. [71 upper panel, for the full statistics of central and mid- 
peripheral collisions (b < 7.5 fm) collected in this experiment. The measured values 
are approximately reproduced by the UrQMD predictions which are significantly 
different for the stiff (7 = 1.5) and soft (7 = 0.5) density dependences. 



6. Elliptic-Flow Ratios 

For the quantitative evaluation, the ratio of the flow parameters of neutrons versus 

101 1 

protons or versus Z = 1 particles has been proposed as a useful observable. 1 ^ Sys- 
tematic effects influencing the collective flows of neutrons and charged particles in 
similar ways should thereby be minimized, on the experimental as well as on the 
theoretical side. These include, e.g., the existing uncertainties of the isoscalar EoS 
which will affect the absolute magnitude of calculated flows and the matching of 
the impact-parameter intervals used in the calculations with the corresponding ex- 
perimental event groups. To demonstrate this kind of insensitivity, the calculations 
were performed with two parameterizations of the momentum dependence of the 
elastic nucleon- nucleon cross section, labelled FP1 and FP2, which dif fer in their 
absolute predictions of V2 at mid-rapidity by w 40% for this reaction F^^ 

In Fig. UJ (lower panels), the results for the ratio with respect to the total hy- 
drogen yield are shown. The calculated ratios exhibit clearly the sensitivity of the 
elliptic flow to the stiffness of the symmetry energy predicted by the UrQMD but 



March 3, 2013 21:38 WSPC/INSTRUCTION FILE elliptic'ijmpc 



12 W. Trautmann and H.H. Wolter 



0.00 



-0.05 

cn 

> 

-0.10 



-0.15 



1.0 

.c 

CN 
> 

\ 

0.5 







1.0 

CN 

> 
\ 

> 0.5 





0.2 0.4 0.6 0.8 1.0 

p t /A (GeV/c) 

Fig. 7. Elliptic flow parameters t>2 for neutrons (dots) and hydrogen isotopes (open triangles, 
top panel) and their ratio (lower panels) for moderately central (b < 7.5 fm) collisions of 197 Au 
+ 197 Au at 400 MeV/nucleon, integrated within the rapidity interval 0.25 < y/y p < 0.75, as a 
function of the transverse momentum per nucleon pt/A. The symbols represent the experimental 
data. The UrQMD predictions for 7 = 1.5 (a-stiff) and 7 = 0.5 (a-soft) obtained with the FP1 
parameterization for neutrons (top panel) and for the ratio (middle panel), and with the FP2 
parameterization for the ratio (bottom panel) are given by the dashed lines (reprinted from Refl^JJ 
Copyright (2011), with permission from Elsevier). 

depend only weakly on the chosen parameterization for the in-medium nucleon- 
nucleon cross section. The experimental ratios, even though associated with large 
errors, scatter within the interval given by the two calculations. Linear interpola- 
tions between the predictions, averaged over 0.3 < pt/A < 1.0 GeV/c, yield very 
similar results 7 = 1.01 ± 0.21 and 7 = 0.98 ± 0.35 for the two parameterizations. 
The error is larger for FP2 because the sensitivity is somewhat smaller (Fig. [7j 
bottom panel). 

This analysis was repeated in various forms. With the squeeze-out ratios v^jv\ 
of neutrons with respect to free protons, similar results were obtained, however with 
larger errors. The study of the impact parameter dependence indicated a slightly 
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Fig. 8. Sensitivity of the neutron-proton elliptic flow difference v^~ p = v% — («2 as defined 
in Eg. PH l to the isospin independent EoS and to the symmetry term in the parameterization of 
Refill as indicated. The widths of the bands represent the uncertainty arising from using either 
a soft (K = 210 MeV with momentum dependence) or a ha rd (K = 380 MeV) version for the 
isospin-independent part of the EoS (reprinted from Ref.G23 Copyright (2011), with permission 
from Elsevier). 



smaller value 7 ~ 0.5 for the mid-peripheral event group, again with larger errors. 
It was also tested to which density region around po the elliptic-flow ratios are 
sensitive, with the result that both, sub- and supra-saturation densities are probed 
with this observable 

In consideration of the apparent systematic and experimental errors, a value 
7 = 0.9 ± 0.4 has been adopted by the authors as best representing the power- 
law exponent of the potential term resulting from the elliptic-flow analysis. It falls 
slightly below the 7 = 1.0 line shown in Fig. [4] but, with the quoted uncertainty, 
stretches over the interval from 7 = 0.5 halfway up to 7 = 1.5. The correspond- 
ing slope parameter is L — 83 ± 26 MeV. The squeeze-out data thus indicate a 
moderately soft to linear behavior of the symmetry energy that is consistent with 
the density dependence deduced from experiments probing nuclear matter near or 
below saturation. Comparing with the many-body theories shown in Fig. [3l the 
elliptic-flow result is in good qualitative agreement with the range spanned by the 
DBHF and variational calculations based on realistic nuclear potentials. 

7. Elliptic-Flow Differences 

In an independent analysis, Cozma has used data from the same experiment and 
investigated the influence of several parameters on the difference between the ellip- 
tic flows of protons and neutrons using the Tubingen version of the QMD trans- 
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port mo del ESI They included the parameterization of the isoscalar EoS, the choice 
of various forms of free or in-medium nucleon-nucleon cross sections, and model 
parameters as, e.g., the widths of the wave packets representing nucleons. The in- 
teraction developed by Das et al. was us ed whi ch contains an explicit momentum 
dependence of the symmetry energy part 142143] j n Fig. [S] the difference of neutron 
and proton elliptic flows for five choices of the x parameter describing the density 
dependence of the symmetry energy is reproduced from this work. It shows the 
large sensitivity to the stiffness of the symmetry energy (different x-parameters) 
and the reduced influence of the isoscalar EoS (widths of the bands) . A super-soft 
behavior of the symmetry energy was confirmed to be excluded by the comparison 
with the experimental flow data. 

Because the flow difference is still proportional to the overall magnitude of the 
elliptic flow predicted by the model, also the ratios of the neutron-versus-proton 
elliptic- flow parameters have been calculated and compared to the FOPI/LAND 
dataPSl The best description is obtained with a density dependence slightly stiffer 
than the x = solution which approximately corresponds to the 7 = 0.5 case shown 
in Fig. SJ It is thus rather close to the UrQMD result 7 = 0.9 which may represent 
an important step towards the model invariance that one would ultimately like to 
achieve, not only because of the invariance with respect to the treatment of the 
nucleon-nucleon cross sections observed in the two studies but also because of the 
explicit momentum dependence of the isovector potentials which is implemented in 
the Tubingen QMD but not in the UrQMD version used. 

8. Conclusion and Outlook 

According to the predictions of transport models, the relative strengths of neutron 
and proton elliptic flows represent an observable sensitive to the symmetry energy 
at densities near and above saturation. By forming ratios or differences of neu- 
tron versus proton or neutron versus hydrogen flows, the influence of isoscalar-type 
parameters of the model descriptions can be minimized. 

The comparison of the results obtained from the FOPI/LAND data for 197 Au + 
197 Au reactions at 400 MeV/nucleon with UrQMD transport calculations favors a 
moderately soft to linear symmetry term with a density dependence of the potential 
term proportional to (p/po) 7 with 7 = 0.9 ± 0.4, compatible with predictions of 
ab-initio calculations. With the Tubingen QMD model, a similar result is obtained. 

The explicit proof that the elliptic-flow ratios are probing the nuclear mean field 
at super-saturation densities is, so far, limited to test calculations made with the 
UrQMD They indicate that the strength of the symmetry energy at densities 
both, below and above saturation, are essential. It will be useful to study this in 
more detail. 

The statistical uncertainty of the existing FOPI/LAND data is larger than the 
investigated systematic effects of model parameters and analysis techniques. More 
significant data can thus be expected from a new experiment delivering a compre- 
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hensive data set with sufficiently high statistical accuracy. Explicitly testing the 
predicted dependences on rapidity and transverse momentum will provide useful 
constraints for the models. Additional observables as, e.g., neutron-to- proton or 3 H- 
to- 3 He yield ratios may serve in resolving ambiguities caused by the effective-mass 
splitting between neutrons and protons. It will be interesting to see first results ap- 
pearing from the recen tly completed measurements of the ASYEOS Collaboration 
at the GSI Laboratory!^ 

Because of the quadratically rising importance of the symmetry energy, the 
continuation of this program with systems of larger asymmetry is very promising 
and important. The lower luminosities to be expected from the use of secondary 
beams and isotopically enriched targets will have to be compensated with efficient 
detector setups. The neutron detector NeuLAND proposed for experiments at FAIR 
will offer a h ighl y improved detection efficiency for neutrons in the energy range 100 

VJA I 

to 400 MeVM-^ 1 This will be essential for extending the program also to reactions 
at lower energies for which significant mean-field effects are predicted for directed 
and elliptic flows 

Stimulating and fruitful discussions with M.D. Cozma, M. Di Toro, W. Reisdorf, 
and with the authors of Ref.^B, P. Russotto, P.Z. Wu, M. Zoric, M. Chartier, 
Y. Leifels, R.C. Lemmon, Q. Li, J. Lukasik, A. Pagano, and P. Pawlowski are 
gratefully acknowledged. 
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